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Abstract 

We calculate the parity-violating amplitudes in the nd interaction with pionless effective field 
theory to LO. Matching the parity violating low energy constants to the DDH coefficients we make 
numerical predictions for parity- violating observables. In particular we give predictions for the spin 
rotation of a neutron on a deuteron target, and target and beam asymmetries in nd scattering. 
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I. Introduction 



Hadronic parity-violation has been traditionally analyzed in terms of potential models; 
specifically the DDH model[l|, which is a parity violating single meson exchange picture 
containing seven phenomenological constants. However, there exist well known discrepancies 
between experimental measurements and the DDH modeli2|. Some of this discrepancy is 
no doubt due to nuclear physics uncertainties, but another source may be the use of the 



model-dependent DDE 
proposed by Zhu et al. 



_3otential. A possible solution to these problems has recently been 
2, y] , restriction of experiments to nuclei with A < 4 so that nuclear 
uncertainties are minimal, and analysis using a model-independent picture via effective field 
theory. At low energies, less than m^/MN ~ 20Me\^, such an approach is provided by 
Pionless EFT (EFTy), which has been extremely successful at low energies in the two- 
body and three-body sector for parity-conserving interactions, including interactions with 

n 

external currentsi4]. At low energies inclusion of parity- violation requires only five additional 
low energy constants (LEC's) in the nucleon-nucleon interaction. These EEC's involve all 
possible isospin structures that mix S and P waves with one derivative and are equivalent 
to the parameters originally posited by Danilovj^. The fact that five LEC's are needed 
has also been specifically shown by Girlanda, who does this by performing a non-relativistic 
reduction of all possible one derivative relativistic parity-violating structures that conserve 
Cp|^. Calculations for parity- violation using EFT methods have heretofore been primarily 



focused in the two-body sector. Such calcu 



scattering and in the process np — )■ d'~^\ 



ations include parity violation in nucleon-nucleon 



Parity-violating EFT calculations have been 



done in the three-body sector using a hybrid approach, wherein the parity- violating potential 
is given by EFTy, but is used with wavefunctions determined by either the variational 
method or by a Fadeev integral equation technique 10- 12|. Such calculations include neutron 



spin rotation and beam asymmetry in nd interactions. Recently a paper by Griesshammer, 
Schindler, and Springer predicted the spin rotation of a neutron on a deuterium target up 
to and including NLO effects in EFTyjlSl]. However, they only included order of magnitude 
estimates for the parity violating coefficients and left open the calculation of other possible 
parity violating observables in nd interactions. In this paper we set out to obtain estimates 
for the parity violating coefficients by matching them to the DDH "best" value estimates. 
As well as calculating the neutron spin rotation on a deuteron target at LO in EFTy, 



we also calculate the beam and target asymmetry at LO in nd scattering. Estimation 
of these observables will then allow one to assess the feasibility of nd interactions as a 
realistic experimental probe for the five PV LEC's. Below we calculate the LO amplitudes 
for S-P mixing in nd scattering due to the two-body parity-violating Lagrangian. (Since 
Griesshammer and Schindler showed that no three-b ody parity-violating force occurs up 
to and including NLO, only five LEC's exist at L0[14].) Predictions are made for PV 
observables and numerical estimates are given based on DDH "best" value estimates. In 
a future publication we shall present higher order corrections. The paper is organized as 
follows. In section II we give the form of the two-body parity-violating interaction. Then in 
section III we show what diagrams are needed at LO and how to calculate them. Section 
IV shows how estimates for parity-violating LEC's can be obtained, and in Section V we 
show how to relate our amplitudes to observables. Finally in Section VI we summarize the 
results. 



II. Two-Body Parity- Violating Interaction 

The leading order two-body parity-conserving Lagrangian in the auxiliary field formalism 
is given by 
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t]N^PiN + h.c. 
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where ti (s^) is the deuteron (singlet) auxiliary field|l5|. Here Pi = -^o'20'iT2 projects out the 
^Si channel and Pa = -^cr2T2'^a projects out the ^5*0 channel. The auxiliary field formalism is 
equivalent to the partial wave formalism in which only nucleon fields are used as can be seen 
by integrating over the auxiliary fields and using a field redefinition {3]. In this formulation 
the two-body parity-violating Lagrangian amplitude is given by a form including five low 
energy constantsis]. 



FIG. 1. Infinite sum of nucleon bubbles contributes to LO deuteron propagator 
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where a V & 



i(^)b—(^a)b, andX = diag[l, 1, —2] projects out the isotensor contribution. 
The deuteron kinetic energy and the term A|q^^^ are sub-leading with respect to '^[^ly 
Thus at LO the bare deuteron propagator is given by i/A^^^^^ which is then dressed by an 
infinite number of nucleon bubbles as seen inf Fig. [1] The resulting propagator depends 
on A(.'^\^ and yj, with values adjusted such that the deuteron propagator has its pole at 
the correct value. A similar calculation can be carried out for the propagator of the singlet 
auxiliary field. This procedure has been carried out in many papers, the end results for the 
LO deuteron and singlet propagator are listed below [l5|. Also we include the constraints 
imposed on the coefficients A^^'^^-^^, yj, and their ^5*0 counterparts at LO. Note the presence 
of the parameter fi, which is a cutoff imposed by using dimensional regularization with the 



PDS subtraction scheme 



171]. (Here 7^ = 45.7025 MeV is the deuteron binding momentum 



and 7t = 1/at, where at = —23.714 fm is the scattering length in the ^5*0 channel. 
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III. Three-Body Parity Violation LO 

As in the parity-conserving case for three-body interactions one needs to solve an infinite 



sum of diagrams for the parity- violating amplitude at LQ|15l. Il8| . leading to a coupled set 
of integral equations. Numerical solution is necessary, as such integral equations cannot 
be solved analytically. In general we must solve a set of four coupled integral equations. 
However, since parity- violation is so small, Gpni^. ~ 10^^ we can ignore second order PV 



terms. Then the integral equations for the parity-conserving amplitudes decouple 191] . and 

n n " 



are exactly the same as in previous papers [15|, ll8|, |20|]. The remaining coupled parity- 
violating integral equations at LO are shown in Fig. [2l where the boxes represent parity- 
violating vertices, the double line the dressed deuteron propagator, the double dashed fines 
tfie dressed singlet propagator, and tfie line witfi arrow tfie nucleon propagator. Tfie tfiick 
lines represent a sum over botfi deuteron and singlet propagators. Tfius tfie tfiick line allows 
one to represent two Feynman diagrams witfi a single diagram. Tfiere are also diagrams 
wfiere two dibaryon lines and two nucleon lines meet at a single vertex, due to tfie tfiree- 
body force term in tfie Lagrangian. 

Tfiis tfiree-body force term enters at LO, only in tfie Doublet S-wave cfiannel. (Note we 
fiave not yet projected out any specific cfiannel.) Tfie momentum integrals are regulated 
using a sfiarp cutoff A. Tfie tfiree-body force term is cutoff dependent. Tfiis cutoff is 



convenient because it can be implemented straigfitforward 



force term in tfie Doublet S wave cfiannel is given by 2l|, |22 1 



y numerically. Tfie tfiree-body 



«(£,A) = ^ (6) 
Tfie cutoff dependence of tfiree-body force term ifn(A) is cfiosen so tfiat tfie Doublet S wave 



amplitude produces tfie correct scattering lengtfi[l8l. l2d . l23l |. (As noted earlier tfiere is no 
need to include a parity-violating tfiree-body force, as it fias been explicitly sfiown tfiat no 
sucfi term exists up to and including NLojl^.) 

To first order in parity-violation tfie integral equation for tfie parity-violating amplitude 
is given by tfie sum of diagrams sfiown in Fig. [31 Tfie first diagram corresponds to a pure 
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PV — \ + PC \ + PV 




FIG. 2. Integral Equations for Parity Violation at LO (Note diagrams where lower vertices are 
parity- violating are not included) 
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FIG. 3. Parity- Violating Diagrams at LO (Note diagrams where lower vertices are parity- violating 
are not included) 

parity-violating transition with no scattering in the initial or final channel. The next set of 
diagrams has a parity- violating transition with scattering either in the initial or final channel 
but not both. (Also note that for these diagrams the singlet field acts as an intermediate 
state, which can only exist in the Doublet channel.) Finally we have the set of diagrams 



with a parity- violating transition and scattering in both the final and initial channels. 

Summing all of these figures one finds the parity- violating amplitude given in Eq. ([7]), 
where k is the incoming nucleon momentum and p is the outgoing nucleon momentum in 
the cm. frame. Since our diagrams are on shell we have |k| = |p|, and the total energy 
in the cm. frame is given by E = — j^. The vector index letter i (j), represents 
the initial (final) deuteron auxiliary field polarization. Finally the Greek index a is the 
initial (final) spinor index and a (b) is the initial (final) isospinor index. 
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The vector v, p.wcts out the uudeou-deuteton amphtude i„ cluster-configuration space 
and is defined as [20] 




and the parity-conserving amplitudes t are a vector defined as follows 



7 
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where t^d^^d is the amphtude for nd scattering and tjq^^m is the amphtude for nd going to 
a nucleon and a singlet combination of the remaining nucleons. (Note that we have not yet 
projected out Quartet or Doublet channels.) The expressions D(£^,q) and (i/C"' )^^(q, ^, go) 
are both matrices defined via. 



, iDd(E,q) 
.D(E,q)= '^^ (10) 

iA(^,q) 



(zK;^'^)ft(q,^,go) ^ ^ X (11) 

iq2 + q-€+i^ +lk^ + ^j-Mj,qo-ie 



X 



(/ci>V^)f„(q,^) (/c|.VOfl(q,^) 



where the functions (^py"")aa (q''^)' which contain all of the parity- violating dependence 
are defined as 



KVfl A P) = y<i9'''-''''{'^W^ + m + W^-'^^e^^^(aV^)^(r3)^(k + 2pY (12a) 

KVXa P) = + 2p)^ + ^y.5(lvJ5°e^^^(<xV^•)^(r^)^(k + 2p)^ 

(12b) 

+ yt/^^-'^Mf (r^)^(2k + p)^- - ty,g'^^-'^^e^''^{T^Ts)l{a'mk + p)^ 
{JClYtl (k, P) = (k + 2py + zy,/^-^^^6'^na'=)^(r3r^)^(k + 2p)^ 

(12c) 

+ y.4v:5°(aV0f (r^)^(2k + p)^ - i^io-:j;°e3^^(aV0f (r^)^(2k + p)^ 

{iCfy^X'a (k, P) = l/*^(K5°(-0^Jr^T")^(k + 2pY + ^y4'Z-:^e^^^i-ra{r^rya{^ + 2p)^ 

(12d) 

+ (r^r^)^(<T0^(2k + p)^ - iy,,;iv:?^^^^(r^r^)^('-0^(2k + p)^ 



(Note that the capital letters A,B, and C are used for the singlet auxiliary field polarization 
and the lowercase letters i,j, and k are used for the deuteron auxiliary field polarization.) 
Integrating over the energy and picking up the poles from the nucleon propagators in our 
diagrams Eq. ([7]) becomes. 



('?.):Ap) = ^vJ(;c'-):(k.p)v, (13) 

where 



(/C^')^a(q,^) = X (14) 

Now that we have derived the parity- violating amplitude, we note that it contains the 

related scattering amplitudes from the parity-conserving sector. Such parity-conserving 

i 

scattering amplitudes are calculated in[20], by numerically solving Fadeev's equation in an 
angular momentum basis. However, as part of this solution one runs into singularities along 
the real axis. To overcome this difficulty the method of Hetherington and Schick is employed, 
in which the axis of integration is rotated into the complex plane, therefore avoiding the 



singularities 24j-|26|. One can then use the solutions along the deformed contour to solve for 
the amplitudes along the real axis. Details of the procedure to calculate these amplitudes 
can be found inj27|. In order to use the solutions to Fadeev's equations we need to project 
out our parity-violating amplitude into an angular momentum basis. However, unlike the 
parity- conserving sector, the PV amplitudes mix different angular momentum states. Also, 
since at leading order, spin and orbital angular momentum mix, the appropriate angular 



momentum basis to use is the total angular momentum J = L + S. Thus we express our 
parity- violating amplitude as 

oo M=J J+S J+S' 
J=0 M=-J L=\J-S\ L'=\J-S'\ S,S' 

where the spin angle functions are given by 

= E cT,'sT'''Yr^i^)xT (16) 

mL,ms 

being the spinor part of the spin-angle functions, C*™^'^'^'^ the appropriate Clebsch- 
Gordan coefficient, and the appropriate Y^^(k) spherical harmonic. Since the spin-angle 
functions are orthogonal, we can project out the amplitudes in our angular momentum basis 
via 

ti^S',Lsik^P) = i^ld^'^l {%%m*tpv{K P)^.55(k) (17) 
At sufficiently low energies S-P mixing will dominate. Thus we will calculate only the 
amplitudes tfgf^g (Note J=S here since L=0) for all possible values of 5" and 5". All spin 
and angle dependence is contained within the matrix (/C*-')f^(q, -6), and the appropriate 
projections in J, L,and S can be found in the appendix. Going to a partial wave basis we 
finally obtain an expression for the parity-violating partial wave amplitudes. 



tpvifs,,i,s(k,p) = ^wjK{k,p)is'.LS^,+ (18) 
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This expression contains the parity-conserving amplitudes in the partial wave basis of total 

— * — * — * 

angular momentum J = L-l-S. (These are equivalent to the parity-conserving amplitudes in 
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the partial wave basis of orbital angular momentum.) It can be shown straightforwardly that 
the parity-conserving amplitudes are independent of total angular momentum J. Thus we 



can use the parity-conserving amplitudes as calculated numerically by|l5l.ll8l.l20l] and perform 
the integration numerically in order to obtain the associated parity-violating amplitudes. 

Before integrating Eq. fITS]) we multiply by the LO deuteron renormalization Zj^ = 
(87r"/w) / (My?/^) [isl. and use the renormalized parity-conserving amplitudes. We find that 
all the parity-violating LEC's occur in the combinations. 

3Si-lPi 3o-3p, n^So-'iPo ^^So-^Po 

Vd ' Vd ' yt ' vt 

(Note g^^^J2) does not appear as a A/ = 2 transition is not allowed for a first order parity- 
violating transition in nd scattering.) For the sake of convenience we find it useful to make 
the following definitions. 

35i-iPi Sc^.Sp, n'-^O-^Po n^O-^'Po n^^O-^^O 

g^' g^' 9(AI=0) 9^AI=1) ^(A/=2) 

gi = ,g2 = ,g3 = ,^4 = ,^5 = 

Vd Vd Vt Vt Vt 

Since these coefficients are unknown, we will write the parity-violating partial wave ampli- 
tudes as follows, where {tpvjjs' Lsi^iP)Y calculated by setting Qj = 0,j ^ i and gi = 1 
in the parity-violating partial wave amplitude. Thus the parity-violating amplitude can be 
written as. 

4 

tpvi^s',Ls{k,P) = {tpvi^s',Lsik,p)y (19) 

1=1 

IV. Parity- Violating Potential 

It is clear from Eq. f|T9|) that in order to obtain numerical values for the parity-violating 
amplitude one needs to know the size of the coefficients gi, which at this time are not 
determined from either theory or experiment. Nevertheless, we can obtain estimates by 
matching the gi onto the familiar DDH coefficients. We will carry out this procedure in 
three steps. First we match the DDH coefficients onto the coefficients of the Zhu potentialP 
Then we match the Zhu potential on to the Girlanda potential 
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12l |. Finally we project 



the coefficients of the Girlanda potential onto the coefficients of the auxiliary field formalism. 
We also show how all these formalisms can be matched to the familiar Danilov parameters 

11 



The DDH model 

mesons vr, p, and u. 



2M, 



N 



(20) 



2M, 



AT 



with the strong couphngs given approximately by g'^j^j^/An ^ 13.5 and gp/4:n = hg^/4: 



.67, while the magnetic moment terms are approximately Xp = ~ = 3.7 and = 
i^p + i^n = —.12. The phenomenological weak interaction Hamiltonian posited by DDH 
consists of seven weak coupling terms 



{21) 



0" ii.uk'^ 



DDH attempted to obtain theoretical predictions for the seven constants using SU(6) sym- 
metry and quark model techniques. However, due to the difficulty of this calculation they 
were only able to come up with reasonable ranges and "best" values as shown in Table [B 
(Also shown are estimates by other groups.) 



Coupling 


DDH[1] 

Reasonable Range 


DDH[1] 
"Best" Value 


DZ[29] 


FCDH[30j 




^ 30 


+12 


+3 


+7 




30 ^ -81 


-30 


-22 


-10 


hi 


-1^0 


-.5 


+1 


-1 


hi 


-20 -29 


-25 


-18 


-18 


hi 


15 ^ -27 


-5 


-10 


-13 


hi 


-5 -f -2 


-3 


-6 


-6 



TABLE I. Weak NNM couplings. All numbers are quoted in units of the "sum rule" value 
Sr = 3.8 X 10"^ 



^ Since CP is conserved there are no neutral pseudoscalar mesons tt", 77, or 77' by Barton's theorem 



28| 
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The form of any parity- violating potential can be written as a sum of operators O^j with 
corresponding coefficients c^, where a refers to the specific potential of interest.. 

V- = J2<0lf (22) 

n 

At the lowest energies the component of the operators that contain momentum is of two 
forms. 

{P..,/n(^^,)} (23) 



X 
X 



(n) 

(n) 
ij- 



where Pij = (pi — P2)/2 is the momentum of the nucleon-nucleon system in the cm. frame. 



TABLE II. Parity- violating potential in DDH and Zhu formalism. Tlj = (3r?T? — Tj • tj). (Note 



A?, ~ 'inFTT is the chiral scale 



31 



32l |. where = 92AM eV is the pion decay constant) 
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„DDH 


^Zhu 
'-n 




fZhu 
J n 


0^"^ 
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1 dnNN f 

2V2Mn 


X 


Mr) 


fm{r) 


{n^T,)^{a, + aj).xf^^_ 


2 


Mn P 


tt2C3 

X 


fp{r) 


fm{r) 


{n ■ Tj){ai - aj) ■ X\f^^ 


3 


9p(l+Xp) 1,0 


tt2C'3 

X 


fp{r) 


fm{r) 


{n-T,)ia,xa,).x\f_ 


4 


2Mjv P 




fp{r) 


fm{r) 


(r,+T,r{a,-aj).x\^l 


5 


9p(l+Xp) Ll 
2Mjv P 


X 


fp{r) 


fm{r) 


iT,+T,)%a,xa^).x[^l 


6 


9P ;^2 
2V6Afjv P 


X 


fp{r) 


fm{r) 


Tijia^-aj)- X[f^ 
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2V6Mn P 


Tree's 

X 


fp{r) 


fm{r) 


Tlj{a,xa,).x'^_ 
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Mn ^ 


tt2Ci 

X 


Ur) 


fm{r) 


{a.-a,).xfl^ 
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2Mn ^ 
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(r,+r,)^(<T,-o^,)-^!r+ 
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2Mn ^ 




Ur) 
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fp{r) 
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(r,-T,)^{a, + a,).x\fl 
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2Mn P 





fp{r) 





{T,XT,)^ia,+a,).x[]'l 
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The coefficients, operators, and regulator functions /, 



DDH I 



for the DDH potential and 



Zhu/ 



for the Zhu potential are given in Table [TTl The functions /, 



DDH , 



Tij) are Yukawa 



functions, where the mass corresponds to the appropriate meson j^. However, at the lowest 
energies the functions for the DDH potential can be written as fi{r) = ;^(5^(r), where 
i = 7T,p, or a; [2]. Likewise the functions fm{r) = 6^{r) for the Zhu potential become delta 
functions in the low energy limit. Thus at low energies the DDH potential and the Zhu 
potential can be trivially matched yielding jsl 
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l + Xp^ 4.7 
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(25) 
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(jDDH 



qDDH 
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A? 



2M^m2 



g^hl 2.25 X 10-^ 



A3 

, ]^ bestguess 



2MNml 



1.35 X 10"^ 



A' 



X ^ bestguess ^^^g ^ ^Q_6 



A? 



2MNml 



g^hl '^^^^^ 7.64 X 10- 



DDH 
6 



-.^^ g hl ""'^^ -7.80 X 10-^ 

' 9.19x10-^ 



"4v^M^m2- 



(26) 



As first pointed out by Danilov, one needs five parity- violating terms at the lowest energies 
in the two-body sector [5|], since only S-P mixing is involved. By conservation of angular 
momentum the state ^6*1, can only connect with the states or ^P\. Since is an isosinglet 
there is a unique way to get to the isosinglet state and isotriplet state ^P^. Similarly, the 
state ^S'o can only connect with the state ^P^. However, both ^5*0 and ^Pq are isotriplet states 
so the operator connecting these states can carry A/ = 0, 1, or 2. The existence of five unique 
operators which characterize parity violation at low energy appears to be in contradiction 
with the DDH and Zhu potential, which involve ten different operators. However, at low 
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energies five of tliese operator structures are redundant as sliown by Girlanda6|. In tfiis 
procedure one begins witli all possible one-derivative P violating CP conserving relativistic 
terms. 



01 = 








C»2 = 
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Using Fierz transformations and the equations of motion, there exist six identities 



03=0i d2 + d4=MM{02 + 0^) (28) 

O2 - O4 =20e 02 - 04 = - 2MM - Oe 

O3 + 3d, =2Mjv(Oi + O3) =MM 

reducing the number of unique operators to six. However, in a non-relativistic reduction it 
turns out that two of the operators have equivalent structures leaving five unique operators 
at the lowest energies. The resulting parity-violating Lagrangian in the Girlanda formalism 
is given by 



£Gir _ ^^(ATt^A^ .nH\/N- N^NNH V -^N) - g^tijkN^aiNVjiN^akN) (29) 

- g2ei,k[N^T3a^NVjiN^akN) + N^aiNWjiN^nakN)] 

- Q^XabeijkN^TaCJiNVjiN^nakN) + g^tabz V {N^TaN) ■ N^naN 

(In Eq. ( I29l) a factor of has been absorbed into the coefficients. This notation agrees 
with the notation of Phillips, Schindler, and Springer [rj].) With this parity- violating La- 
grangian one can compute the Girlanda potential which takes on the following form given 
by Eq. (!22|) . where n runs from one to five. 
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TABLE III. Parity violating potential in Girlanda formalism. Tij = (3t?t? — Ti ■ tj). 



3 -2^5 S^{r) Tijiai - aj) ■ Xf.l 

5 2^1 6^{r) {a,xaj)-x\f^_ 



Using ( !28l) one can reduce the Zhu potential to a set of five operators, allowing the matching 
of the Zhu coefficients onto the Girlanda coefficients as shown in Table IVl 

For our calculations, we also require the coefficients in the auxiliary field formalism Eq. 
([2]). This matching of the Qi and Qi requires two steps. One first performs Gaussian integra- 
tion over the auxiliary fields followed by a field redefinition to rewrite the Lagrangian, Eq ([2]) 
in terms of nucleon fields, as done by Schindler, and Springer jsl]. Then one can match this 
partial wave formalism onto the Girlanda formalism by performing Fierz rearrangements 
and using the constraints Eq. ( l28|) . yielding the results in Table IVl (This has also been done 
using a different method by Phillips, Schindler, and Springer j?!.) 

Finally we wish to match the Girlanda potential onto the Danilov potential which is given 
by Eq. ( I22l) in Table \1V\ where n runs from one to five. 

TABLE IV. Parity violating potential in Danilov formalism. Tij = (^Tfr? — Ti ■ tj), Pq = 

i (1 - (Ti • ffj), Pi = i (3 + (Ti • ffj). 



1 iPthd ^SHr) {r,-r,na,+a,).x[]l 

2 h^lht ^SHr) {n+r,r{a.-a,).x[f^ 

4 Ai/7d ^SHr) {B,-a,)P,.xfl^ 

5 AO/7* ^5^{v) {a,-a,)P,.xfl 
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In order to match the Girlanda formahsm to the Danilov formahsm we note the identity 

{P\HV,6'{r)]\S) = {P\ {-zV,6'{r)} \S) (30) 

which follows trivially since P waves are zero at the origin. Next we make use of the identical 
identities in spin and isospin space. (Note PJ" = ^ (1 — ■ t^-) and -PT = ^ (3 + ■ 7^-)) 

z(«T, X a,) = {a., - aj) (Po - ^0 (3i) 

z(r.xr,r = (f.-T,r(Po--Pn (32) 

The isospin operator i(Tj x Tj)^ only appears with the spin operator (<Tj + ct^) in the 
Girlanda potential. Since this spin operator only projects out the triplet state of the S 
wave, only the isosinglet part of the operator i{Ti x tj)^ is projected out. Thus by Eq. 
( 132|) we find that in combination with the spin operator (ctj + <Tj) the isospin operator 
i{Ti X Tj)^ = (tj — Tj)^. Finally using Eqs. ( 130|) . ( 13T]) . and the fact that the identity / 
is / = Po + Pi one can straightforwardly match the Girlanda coefficients to the Danilov 
coefficients, giving the results shown in Table |Vl Also shown in Table |V] are the relation 
between the Zhu, Girlanda, Auxiliary, and Danilov formalisms. The primary goal in low 
energy hadronic parity violation is to determine the value of the Danilov parameters. At low 
energies all of these different EFT formalisms can be shown to be equivalent to the Danilov 
parameters, as shown in Table |Vl Thus one can use whichever formalism is more convenient. 

TABLE V. Translation between various formalisms of parity violating potential 



Zhu 


Girlanda 


Auxiliary 


Danilov 


^ (Ci -C,- 3{C, - (73)) 
^ (Ci + Ci + (C3 + Ca)) 


tog^ 

MN-rtVSg 


-2V2gi 

-4V2g2 

-2^/2^3 

-V2g4 

-4^355 


Xt 
Pt 
A° 
A.^ 
A^ 



17 



Having matched the auxihary coefficients gi to the Zhu coefficients we can now use the 
matching of the Zhu coefficients to the DDH "best" values to obtain estimates for the 
auxihary coefficients which yields. 



91 

92 
93 
94 

95 



4^27? 



9uiXui ,0 
9-kNN 



•^9pXp , 
Mjvm2 P 



4y2M, 



fn- 



9p 



-7.89 X lO^^MeV"^ 
9u 



2 P 



9u;{2 + Xu) ^o ^p(2 + Xp) ^o 



9p{^ + Xp) ,i 9u{'^ + Xw) ,i 



8^2 



TT 



9pi'2 + Xp) j^2 
VQMNml " 



(33a) 

1.43 X 10-^°MeV-^ (33b) 

(33c) 

1.67 X IQ-^^MeV"^ (33d) 
2.50 X IQ-^^MeV"^ (33e) 



4M^m2^^_ 
- 8.53 X IQ-^^MeV^^ 



V. Spin Observables 

Having calculated the various parity-violating amplitudes, we can now relate them to PV 
observables. One such observable is the neutron longitudinal asymmetry Aj^ 33|. In this 
case we scatter longitudinally polarized neutrons from an unpolarized deuteron target, and 
measure the difference of the two cross sections. 



A. = "-±-^ (34) 

Here (T+ (o"_) represents the cross section of positive (negative) helicity neutrons. 

In order to calculate observables, we need to write them in terms of the partial wave 
amplitudes calculated above. We denote the transition matrix by the operator M. Of 
course, M is not diagonal in the orbital angular momentum basis, but rather is diagonal in 
terms of total angular momentum. Defining Mm'^,m,'2;m,i,m2 the T matrix where the nucleon 
has initial (final) spin m2, (^^2), and the deuteron has initial (final) spin mi, (m'^),it can be 
shown 



Mm'„m',.,^,,m, E E E E v/2LTTC™'™^ (35) 

J L,L' S,S' m,s,m'g m'^ 
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Observables are most easily written in terms of this matrix Mjn[,m'2;mi,m2- Having Eq. (135|) . 

,m2;mi,m2 



which gives M^'^ m^;mi,m2 ill terms of the calculated functions Mfigi j^g, we can calculate 



observables by truncating the sum over J, L, and L' at some reasonable level. Thus the 
observable A^r is given by 



A 



m\ ,mU mi, 7712 



N 



(36) 



m2;,m2,m2 I 



Explicitly summing over values of angular momentum from to 1, and spin and J from 
J = 1/2 to J = 3/2 we find. 



Am =-Re 



11/2 



03/2,03/2 ^ ^^^13/2, 13/2 j I ^^^^ 13/2,03/2 



H^-'01/2,0l/2l + ^H'^03/2,03/2l +'^^'^11/2,11/21 



'03/2,03/2 I 



'13/2,13/2 1 



(37) 

/ 



By the optical theorem we note that An can also be written as 



A 



N 



Em (Mn,l/2;m,l/2 


6=0 — ^m,-l/2;m,-l/2 


e=o) 


Em (Mn,l/2;m,l/2 


e=0 + Mfn-l/2;m.,~y2 


9=o) 



-Im 
3 

Im 



11/2,01/2 
^^^01/2,01/2 



2 V2<;2,oi/2 - 4M;f 



-/2,03/2 



2v^M;//;oi/2 



2M, 



+ 3M;//\,,/, + 3M- 



03/2,03/2 



L/2 
'13/2,13/2 



(38) 



Another parity- violating observable is the spin rotation of the neutron as it passes through 
a deuteron target. In this experiment the neutron is transversely polarized and the rate of 
change of the rotation angle with respect to the distance traveled is 10 



12 



33|. 



dcf) 
dz 



AMnN 



9k 



^Re[M, 



m,l/2;m,l/2|9=0 



m, — i/2;m,— 1/2 1 



^0] 



(39) 



where N is the number of scattering centers per units volume, and k is the momentum of 
the neutron in the cm. system. Using Eq. (|35|) we find 
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dz 



AMnN 
27k 



Re 



^^-'11/2,01/2 ^ ^^^^^^13/2,01/2 ^^^-'11/2,03/2 ^VOJW^ay, 



2,03/2 



(40) 



The final parity-violating observable that we will consider is the deuteron target asym- 
metry. In this case an unpolarized beam of neutrons is scattered from a polarized deuteron 
target. At first the deuteron target is polarized in the positive z direction, where z is the 
direction of the neutron's initial momentum. Then the neutron scatters off the deuteron 
target and we measure the cross section ai. In addition the deuteron target is polarized in 
the opposite direction and we measure the cross section ct_i. If there is parity violation one 
will find that o"i 7^ cr_i, and we define the target asymmetry as 



A 



D 



^1 - CT-l 

ai + (7_i 



(41) 



A 



mj/jfrij "^2 



|-'"m'j^,m2; — l,rri2 I / 



D 



,m2;l,m,2 

JTij^/ ,m'2 "^2 



(42) 



- Re 



^ \^^^-'ll/2,lV2 ^ ^^-'01/2,01/2 ) 1^^*^11/2,01/2 ) ^ V Z yVly^i^ -^^i^ -h JWoi/2,01/2 ) 1^^*^13/2,01/2 



-2 f m'/' + m'/' 

1^*^11/2,11/2 + ^'^03/2,03/2 



11/2,03/2 



y2 



13/2,13/2 ^ ^"03/2,03/2 / I ^"13/2,03/2 



3/2 



i<;2,oi/2i' + 2|Mo3/2,o3/2r + m%,,.i.? + mti..x.i.? 



Again Ad can also be calculated via the optical theorem. 



A 



D 



_Em2l"^(^l,'"2;l,m2 


6>=0 — M_i^m2;-l,m2 


e=o) 


Em2 (^I,m2;l,m2 


61=0 + ^-I,m2;-l,m2 


6l=o) 



(43) 



Im 



2m;//\ov2 + ^<;2,oi/2 - '^^Zm. + 2v^^it;2,o3/2^ / 



Im 



^^^01/2,01/2 ^ ^^*^03/2,03/2 



1/2 



VI. Results 

Plotting our results for beam and target asymmetry as a function of center of mass 
momentum k we find the plots given in Fig. HI The thickness of the plot denotes the cutoff 
variation. The cutoff variation runs from 200 MeV to infinity. It appears that the results 



begin to converge after 900 MeV as found by other authors 
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20[ | . Also it should be noted that 



the cutoff variation for the beam and target asymmetries at low energies is actually smaller 
than as shown in the plots, and is displayed with the given thickness in order that the plot 
be visible. The cutoff variation of the spin rotation which is not plotted, is more significant 
than for Ajy and A^, which should come as no surprise since we are dividing by a factor of 
momentum in Eq. f l40p . (We use a liquid deuterium density of = .4 x 10^^ atoms cm~^ 
This means that the spin rotation is sensitive to higher order effects and, since we only did 
a LO calculation, the cutoff sensitivity should be much greater. The plots for the beam and 
target asymmetries extend all the way to 50 MeV. However, the observables cannot be taken 
seriously at these high energies as higher partial waves and higher order contributions will 
become important. (It should also be noted that a significant difference was found at about 
20 MeV if P waves were not included in the parity-conserving amplitudes.) The spin rotation 
observable at zero energy ranged from 6.57 x 10~^ rad cm"^ to 6.82 x 10~^ rad cm~^ due to 
cutoff dependence. This order of magnitude falls in line with previous calculations 10l-ll2|. 



Also the beam asymmetry ranges from 6.79 x 10^^ to 6.83 x 10^^ due to cutoff dependence, 
at Eiab = 15 keV, and this also agrees well with previous calculations [11]. Finally the target 
asymmetry varies little due to cutoff dependence at Eiab = 15 keV with a value of roughly 
7.84 X 10"^. 



Beam Asymmetry 



Target Asymmetry 





FIG. 4. Beam and target asymmetries as function of cm. momentum k 



Finally in order to compare with possible experiments a table of all three observables in 
terms of their contributions from each of the Qi is given below. The spin rotation is given 
at zero energy and the beam and target asymmetry energy of 15 keV. 

In order to obtain a prediction for the observable each row is multiplied by the appropriate 
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value of Qi and then these products are added together to yield the observable. Below are 
two tables with different values for the cutoff. The first table shows the cutoff at 200 MeV 
and the second at setting the cutoff to infinity. 



TABLE VI. value for cutoff A = 200 MeV 



9i 


Rotation, Eiab = 


keV 


An, Eiab = 15 keV 


Ad, Eiab = 15 keV 


1 


-9.33 rad cm-^ 


MeV 


-14.4 MeV 


8.92 MeV 


2 


-18.1 rad cm~^ 


MeV 


-39.6 MeV 


-59.7 MeV 


3 


-5.11 rad cm~^ 


MeV 


-1.83 MeV 


1.65 MeV 


4 


3.40 rad cm~^ 


MeV 


1.22 MeV 


-1.10 MeV 



TABLE VIL value for cutoff A = oo MeV 



9i 


Rotation, Ei^b 




keV 


An, Eiab = 15 keV 


Ad, Eiab = 15 keV 


1 


-9.58 rad cm" 


-1 


MeV 


-14.5 MeV 


9.15 MeV 


2 


-19.0 rad cm" 


-1 


MeV 


-39.9 MeV 


-59.8 MeV 


3 


-8.35 rad cm" 


-1 


MeV 


-2.71 MeV 


2.47 MeV 


4 


5.56 rad cm~ 


1 


MeV 


1.81 MeV 


-1.65 MeV 



Finally we should note that the calculation for the beam and target asymmetries were 
done using both the standard cross section methods and the optical theorem. Plotting the 
results from both, we found they were indistinguishable. This agreement confirms that our 
amplitudes are unitary and acts as a check on the validity of our results. For the values 
quoted in the table, it was found that for the beam and target asymmetries, as well as the 
values from either the cross sections or optical theorem agreed to less than one percent. 

VII. Conclusion 

Above we calculated the low energy parity-violating nd transition amplitudes using pion- 
less EFT. Matching the auxiliary field formalism onto the DDH potential, we made predic- 
tions for the coefficients of the auxiliary field formalism. Using these amplitudes, we were 
able to make predictions for the spin rotation, beam asymmetry, and target asymmetry in 
low energy nd interactions. The values obtained for the neutron spin rotation and beam 
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asymmetry are in agreement with values found by other authors 10l-ll2j . Unfortunately due 
to the smallness of these values they will require very precise experiments to measure. How- 
ever, the five LEC's used in parity-violation at LO are not very well determined. Thus 
upon further experiments it may be found that the values of the LEC's are such that the 
observables for nd interactions are larger than predicted. 

The largest contribution to parity violation was shown to come from the coefficient ^2, 
which contains the pion contribution, and such experiments should then allow one to deter- 
mine its value. It is noted that to first order in parity violation the A/ = 2 ((75) term does 
not contribute. Thus nd scattering is sensitive to four out of the five PV coefficients. 

In principle we should be able to calculate to NLO in EFTy without the need for parity- 
violating three-body forces[l^. Griesshammer, Schindler, and Springer calculated the NLO 
parity violating amplitudes using the partially resummed approach which introduces higher 
order contributions at NL0|13]. However, to calculate the NLO contributions without higher 
order contributions, one must calculate the full off shell LO amplitude. Since a calculation of 
the full off shell LO amplitude is numerically expensive it will be left to a future publication. 

A. Appendix 

In projecting out the amplitudes one has to perform angular integrations which are given 

by 



UjL =^ J dn, J dfi^_l_F;^^*(p)y--(k) (A;Fr(k) + 2pFr(p)) = (A.l) 



a + k ■ p 

2iv -|- 1 Q Q ^m,L,m,m'y 



where 



Ql = U' ^dx (A.2) 
2 J_i X + a 

are functions related to the Legendre polynomials of the second kind up to a factor of (— 1)'^, 
and Pl{x) are the standard Legendre polynomials. 
The projection we must carry out is of the form 

m,P)i's',Ls = ld^>^l Wl'5'(p))* (/C^Ofa (k,p)3^fL5(k) (A.3) 
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(Note the polarization and spin indices are summed over corresponding indices that are not 
exphcitly shown in the spin angle functions). Each matrix element of (k, p) has a 

different projection. Each of these four different projections has in turn four pieces given by 
Eqs. f[T^ . Fortunately many of these terms can be related by time reversal simplifying the 
projection considerably. For simplicity we only show how to project out the first g^^^^'^^^ 
piece of the matrix element ('C-'*)^^ (k, p) as given in (fT^ and fll2ap . and simply quote 
the other results. In order to project this term out we use properties of spherical tensors and 
the Wigner-Eckart theorem to reduce the following expression to a sum over Clebsch-Gordan 
coefficients. 



WjL =26^^^(1/2, m'2|aV-'|V2, m2)(k + 2p)^ (A.4) 

= vl 5Z E ^'2^^riT''"'^itr^ 

m,m' k,q 

Using the above expression with fll4p . fll2ap . and (IA.3P and, for the time being ignoring the 
isospin, we find the expression for our projection is. 

^IM2,S' '-'1,1/2, 5 ^L,S,J ^L',S',J 

Y;:'-\p)Y^-{k) (A:rf-(k) + 2prf-(p)) (-1)- 

vl E E v/2^^n,f '"'^ilf'"^'^;:f;'^(V2iiT.iiv2) 

m,m' k,q 

Integration over the angular variable can be carried out trivially by using ( lA.ip leaving a 
sum of products of Clebsch-Gordan coefficients . Then using symmetry properties of the 
Clebsch-Gordan coefficients we find. 



X 
k 



Y,y'ss'L'k{y2\\Tk\\y2) 



1/2^1/211111 

1 S' 1 I 1 1/2 s 
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in terms of 6-j symbols 3^ . 

The sum over k can be removed by use of the identity 351 ] . 



1/2 k 1/2 
1 S' 1 




(A.7) 



yielding 



f 1 1 1 1 f L' 1 L 1 1 
[ 1/2 5 J [ S J J 

Using time reversal invariance and including our coefficients, as well as the isospin pro- 
jection, we find the projection for the ['^(A;,^);^,^,^^^],^^ term is 



' L J L' j 

(A.9) 

" ^ L' J L \ kP 



1 1 1 M L' 1 L 
^2 S' S \ \ S J S' 



1 1 1 M L 1 L' 
^2 S S' \ \ S' J s 



X 



X 



X ^{2kQL'{a)+pQL{a)) 

Now combining isospin with our spin projections we find the \j^{k,p)'l,g,j^g\^^ term pro- 
jected out is 
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L' 1 L \ 1 



\ S J S' \ 

Now using time reversal symmetry we see the projection of the [J^{k,p)'[,g, ];^s\^^ term is 

" [ 5' J 5 J 

+ yi/^-'^^247r4=(-l)''-'''-^5s'V.^C^riL'| ^ ^ ^ { \ ^ 1 L' I 1 ^^^^^^^^ ^ 2pg4a)) 



^ ^ ^'^'^ [1/2 3^2 j[S' J S j kp' 

Now combining the spin with the isospin projections we find the projection of the 



(A.12) 
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